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Abstract 



We study the combined effect of cubic anisotropy and quenched uncorrelated 
impurities on multicomponent spin models. For this purpose, we consider the 
field-theoretical approach based on the Ginzburg-Landau- Wilson y? 4 Hamil- 
tonian with cubic-symmetric quartic interactions and quenched randomness 
coupled to the local energy density. We compute the renormalization-group 
functions to six loops in the fixed-dimension (d = 3) perturbative scheme. 
The analysis of such high-order series provides an accurate description of the 
renormalization-group flow. The results are also used to determine the critical 
behavior of three-dimensional antiferromagnetic three- and four-state Potts 
models in the presence of quenched impurities. 

PACS Numbers: 75.10.Nr, 05.70.Jk, 64.60.Ak, 75.40.-s 
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I. INTRODUCTION AND MAIN RESULTS. 



The critical behavior of systems with quenched disorder is of considerable theoretical 
and experimental interest. A typical example is obtained by mixing an (anti)-ferromagnetic 
material with a nonmagnetic one, obtaining the so-called randomly dilute magnets. They 
are usually described by using the Heisenberg Hamiltonian with a random-exchange term 

HH,r = - J 12 PxPv ^ ' ^> 

<xy> 

where s Xj j are M-component spins and p x are uncorrelated random variables, which are 
equal to one with probability p (the spin concentration) and zero with probability 1—p (the 
impurity concentration). The pure system corresponds to p — 1. The critical behavior of 
these systems is well established, both theoretically and experimentally, see, e.g., Refs. [1-5] 
and references therein. In particular, a new random Ising universality class describes the 
critical behavior of the random Ising model (RIM) above the percolation threshold of the 
magnetic atoms. 

The 0(M)-symmetric Hamiltonian (1.1) is a rather simplified model for real magnets. In 
particular, it does not take into account the presence of nonrotationally invariant interactions 
that have only the reduced symmetry of the lattice and that are due, e.g., to the spin-lattice 
coupling and to dipole-dipole interactions. In this case, a more realistic Hamiltonian is 

H u , c = -J & ■ s y + s*i> ( L2 ) 

<xy> x i 

where a is the anisotropy coupling. For many materials a/ J is relatively small and thus one 
usually neglects these additional interactions. However, this is fully justified only if they 
are irrelevant in the renormalization-group (RG) sense. This issue may be investigated by 
considering the cubic-symmetric <p A Hamiltonian [2] 

r ( 1 1 1 r i2 1 M 1 

H c = j /^|-(^0( a ;)) 2 + -r0(a;) 2 + -t;[0( a ;) 2 ] +^E^) 4 |, (1.3) 

where <f>(x) is an M-component field and r oc T — T c . Analyses of high-order perturbative 
expansions (see, e.g., Refs. [6-10]) show that the cubic-symmetric quartic w-interaction is 
relevant for M > 3, and in particular, for the physically relevant case M — 3. In this case, 
the nature of the transition depends on the sign of the coupling w: if w > 0, the critical 
behavior is described by a new fixed point with reduced cubic symmetry, while, for w < 0, 
the RG flow runs away to infinity, and the corresponding system is expected to undergo 
a weak first-order transition. In the two-component case, the 0(2)-symmetric fixed point 
is stable with respect to the w-perturbation, and thus, if the transition is continuous, it 
belongs to the XY universality class. 

It is of interest to study the effect of quenched disorder on cubic magnets. As discussed 
in Ref. [12], the critical behavior of these materials should be described by the effective 
Hamiltonian 

n H ,rc = - J PxPv Sx-Sji + oEE S i,i + 1212 D x,ijSx,iS x ,j, (1.4) 

<xy> x i x ij 
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where, beside the random-exchange term, a random-anisotropy term is present. Here the 
anisotropy term D xi j is a random quantity that is traceless and has zero average. Note 
that, for small anisotropy and weak disorder, this additional term should be smaller than 
the other ones, being, loosely speaking, proportional to the product of (1 — p) and a, i.e., 
it is a second-order perturbation of the Heisenberg Hamiltonian. This argument is not 
fully justified at criticality, since, if random anisotropy is relevant, it will eventually change 
the critical behavior. Nonetheless, we expect — as we shall see, experiments confirm this 
assumption — a large preasymptotic region in which such term can be neglected. For this 
reason in this paper we will not consider the random-anisotropy term and we shall discuss 
the critical behavior of the model with Hamiltonian (1.4) with D xi j = 0. 

If only the random-exchange term is present, i.e. D x i j = 0, the critical behavior of the 
model (1.4) can by studied using the field-theoretical Hamiltonian 



1,„ 1 XO 1 x ,1 1 



H rc = I d d x\ -{d^{x)f + -r<P{xf + -^{x)<P{xf + -v 0(x) 2 + -w £ 0,(x) 4 \ , (1.5) 



2 1 



a; 



4' 

^- i=i 



where ip(x) is a spatially uncorrelated random field with Gaussian distribution coupled to 
the local energy density. Using the standard replica trick, one obtains the Hamiltonian 
[13,14] 

H e = jd d xCZ\ [(d^ a ,i) 2 + <J + E I (« + + <Aj } > (1-6) 

[ i,a ij,ab ' J 

where a,b — 1, ...M and i,j = 1,...N. The original system, i.e. the randomly dilute M- 
component cubic model, is recovered in the N — > limit. [15] The coupling u is negative, 
being proportional to minus the variance of the quenched disorder. 

The study of the effective Hamiltonian 7i e in the limit N — > provides also information 
on the critical behavior of the randomly dilute antiferromagnetic g-state Potts model for 
q = 2 and 3, with Hamiltonian 

Hdq = J E PxPy$s x , Sy , (1.7) 

<xy> 

where J > 0, s x — 1, . . . , q, and p x are uncorrelated random variables. Indeed, as argued in 
Refs. [20,21] using RG arguments, the critical behavior of the antiferromagnetic three- and 
four-state Potts models on a cubic lattice should be described by the cubic Hamiltonian 7i c 
with M = 2 and M = 3 respectively and with w < 0. The same correspondence holds in 
the random case. The randomly dilute three- and four-state models are respectively related 
to the two- and three-component model with Hamiltonian (1.6) in the limit N — ■> and for 
w < 0. 

Since disorder is coupled to the local energy density, one can use the Harris criterion 
[11] to predict the critical behavior of the model. It states that the addition of impurities 
to a system that undergoes a second-order phase transition does not change the critical 
behavior if the specific-heat critical exponent a P ure of the pure system is negative. If a pU re is 
positive, the transition is altered. This occurs in the Ising case (M = 1), where the addition 
of impurities leads to a new random Ising universality class (RIM). In pure M-component 
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cubic models (M > 1) the specific-heat exponent a pur e is negative; therefore, according to 
the Harris criterion, the pure fixed point is stable against disorder. Nonetheless, disorder 
may still have physical consequences. For instance, it may give rise to new fixed points 
or change the attraction domain of the pure stable fixed point. Systems that are outside 
the attraction domain of the fixed point in the absence of disorder, and therefore show 
a fluctuation-driven first-order transition, may undergo a second-order transition in the 
presence of disorder. Such phenomenon, usually referred to as softening, is well-understood 
in two-dimensional random-exchange models in which disorder is coupled to the local energy 
density. Indeed, it was argued in Ref. [16], and later put on a rigorous basis [17,18], that in 
two dimensions thermal first-order transitions become continuous in the presence of quenched 
disorder coupled to the local energy density. For the cubic model (1.3) for v < 0, it was 
shown [19] that such a softening persists in 2 + e dimensions, while it is absent near four 
dimensions, see the e-expansion analysis of Sec. III. Thus, it is interesting to address this 
issue in three dimensions, where the analysis of Ref. [16] shows that the occurrence of 
softening may depend on nonuniversal features of the model. 

In order to study the RG flow of the effective Hamiltonian Ti, e , we consider the fixed- 
dimension perturbative method in d = 3 and compute the RG functions perturbatively to six 
loops. The analysis of such series allows us to determine the RG flow. We briefly anticipate 
the main results of our analysis. The stability of the stable fixed points of the pure theory 
predicted by the Harris criterion is confirmed. The region v < for any M and the region 
w < for M > 3 is outside the attraction domain of the stable fixed point for all physical 
values u < 0. Moreover, for any M there exists a fixed point in the RIM universality class, 
which is weakly unstable (i.e. with a very small crossover exponent) and may give rise to 
observable crossover effects in physical systems. We do not find fixed points in the region 
v < for any M and in the region w < for M > 3. Therefore, no softening is expected, 
at least for sufficiently low impurity concentration to justify the field-theoretical approach. 
As for the three-state random Potts model, it is expected to have an XY transition as in 
the pure case [22] or a first-order transition, depending on the value of the effective negative 
coupling w. The four-state random antiferromagnetic Potts model is expected to undergo a 
weak first-order transition. 

We predict that cubic magnets with small positive anisotropy have a critical behavior 
controlled by the pure cubic fixed point, which has critical exponents very close to the 
Heisenberg ones. Therefore, experiments should effectively observe the standard 0(3) critical 
exponents. This is in good agreement with the experiments that observe in most of the cases 
0(3) behavior with good accuracy, see, e.g., Refs. [4,23]. On the other hand, systems that 
tend to magnetize along the cubic axes, should show a first-order transition, as in the pure 
case. Experimentally, such a transition has never been observed, probably because of the 
smallness of the cubic anisotropy: due to the very small crossover exponent, ~ 0.01, see, 
e.g., Refs. [7,4], the cubic breaking can be observed only very close to the critical point, i.e. 
for |T — T c \/T c <C 10~ 4 , which is the limit of most of the experiments. As already observed 
in Ref. [24], a posteriori, the experimental results also confirm the validy of neglecting the 
random-anisotropy term in Eq. (1.4) in the experimentally relevant range of parameters. 
Indeed, its presence would give rise to a crossover to a first-order transition in all cases 
[12]. The associated crossover exponent is = 2<f) H — 2 + an ~ 0.39 where 4>h and an 
are the quadratic-anisotropy and specific-heat exponents for the Heisenberg model. Such an 
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exponent is sizable and thus one should have been able to observe the random anisotropy if 
it were not very small. 

The paper is organized as follows. In Sec. II we discuss some general properties of the 
RG flow in three dimensions. Sec. Ill analyzes the RG flow near four dimensions. In Sec. IV 
we present the computation and the analysis of the fixed-dimension pertubative expansion to 
six loops. In the Appendix we compute the differences between the three-component cubic 
and Heisenberg critical exponents, by a reanalysis of the fixed-dimension six-loop expansion 
of Ref. [7] for the cubic Hamiltonian He- 



ll. GENERAL CONSIDERATIONS ON THE RG FLOW 

In this section we discuss some properties of the RG flow of the Hamiltonian (1.6) for 
N — > 0, using general arguments and known results holding for the special cases in which 
one of the quartic couplings vanishes. 

The RG flow in the plane u = is that of the cubic Hamiltonian (1.3) (see, e.g., Ref. [4] 
for a recent review of results). Indeed, for u = the Hamiltonian (1.6) describes N decoupled 
M-component cubic-symmetric models. Therefore, in the plane u = there are four fixed 
points [2]: the trivial Gaussian one, the Ising one in which the M components of the field 
decouple, the 0(M)-symmetric and the cubic fixed point. The Gaussian fixed point is always 
unstable, and so is the Ising fixed point for any number of components M [2,25]. On the 
other hand, the stability properties of the 0(M)-symmetric and of the cubic fixed point 
depend on M. For sufficiently small values of M, M < M c , the 0(M)-symmetric fixed 
point is stable and the cubic one is unstable. For M > M c , the opposite is true: the RG 
flow is driven towards the cubic fixed point. Figure 1 sketches the flow diagram in the two 
cases M < M c and M > M c . High-order perturbative computations in the e-expansion and 
in the fixed-dimension field-theoretical frameworks show that 2 < M c < 3; more precisely, 
M c PS 2.9 [6-8]. This means that the critical behavior of the two-component cubic model is 
described by the 0(2)-symmetric fixed point and therefore belongs to the XY universality 
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class. If M > M c , the cubic anisotropy is relevant and therefore the critical behavior of the 
system is not described by the Heisenberg isotropic Hamiltonian. 

In the three-component case the cubic critical exponents u c , rj c differ very little from those 
of the Heisenberg universality class. Indeed, the analysis of the six-loop fixed- dimension 
expansions of Ref. [7] reported in the Appendix provides the following estimates for their 
differences 

v*-vh = -0.0003(3), Vc~Vh = -0.0001(1), 7c - lH = -0.0005(7). (2.1) 

Note that these differences are much smaller than the typical experimental errors, see, e.g., 
Ref. [4] for a list of experimental results, so that, at present, cubic effects are experimentally 
negligible. Using the accurate estimates of Ref. [26] for the Heisenberg exponents and 
Eq. (2.1), one obtains 

v c = 0.7109(6), r] c = 0.0374(5), 7 C = 1.3955(12), (2.2) 

which are consistent with, but much more precise than, the results obtained from a direct 
analysis, see, e.g., Refs. [7,6]. 

The stability of the pure fixed points against the w-perturbation can be inferred by using 
general arguments [11,25,2]. Since the -u-interaction is the sum of the products of the energy 
operators of the different cubic M-component models, the crossover exponent associated 
with u is given by the specific-heat critical exponent a independently of N, and thus also 
for N — > 0. Therefore, the pure stable fixed point is stable with respect to random dilution 
for any M > 2, since the specific-heat exponent is always negative. For example, for M — 2, 
where the stable fixed point is the 0(2)-symmetric one, we have [27] «xy = —0.0146(8); for 
M — 3, where the stable fixed point is the cubic one, a c = —0.133(2) using Eq. (2.2). 

For v — the Hamiltonian (1.6) describes an MiV-component model with cubic 
anisotropy. The RG flow for N — > is shown in Fig. 2. It is characterized by the presence 
of two stable fixed points. The one for u > 0, w = is in the self-avoiding walk (SAW) 
universality class, but it is irrelevant for our problem, since it is unreachable from the phys- 
ical region u < 0. The other one belongs to the region u < 0, w > and it is in the RIM 
universality class. See, e.g., Refs. [28-32] for recent studies of the critical properties of RIM. 

In the case w — 0, the Hamiltonian (1.6) describes N coupled M-vector models, and it 
is also called MN model [2]. See, e.g., Ref. [4] for a recent review of results. The RG flow 
for M > 2 and iV — > is shown in Fig. 2. Again, the flow is characterized by two stable 
fixed points: the SAW and the 0(M)-symmetric ones. 

For M = 2 and generic N the Hamiltonian (1.6) is invariant under the transformation 

(0M » 02,*) > ^^ 1,< + ^ ' ^ M ~ ^ 2,< )' 

(u , v , w ) — > (u , v + §u>o, -w ). (2.3) 

For N = this transformation maps the Ising fixed point into the cubic one, and the RIM 
fixed point into a new one belonging to the region with -u<0,f>0,w<0. Of course, 
corresponding fixed points describe the same critical behavior. 

In conclusion, the above-reported considerations show the presence of at least seven 
fixed points for M > 3 and eight for M = 2 using the above-mentioned symmetry. Of 
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course, other fixed points may lie outside the planes u — 0, v — 0, w — 0. Thus, a more 
general analysis for generic values of the quartic couplings is necessary in order to obtain a 
satisfactory knowledge of the RG flow. In particular, this would allow us to investigate if 
random dilution may cause a softening of the first-order transition predicted for pure systems 
that are outside the attraction domain of the stable fixed point. For example, systems with 
v < might have a continuous transition if there exists a fixed point in the region u < 0, 
v < 0, while, for systems with M > 3 and w < 0, softening to a continuous transition 
requires a fixed point in the region u < 0, w < 0. 



III. THE RG FLOW NEAR FOUR DIMENSIONS 

The RG flow of the randomly dilute cubic-symmetric model can be investigated near four 
dimensions using the first nontrivial terms of the expansion in powers of e = 4 — d. Using 
the results reported in Refs. [13,33,9], one may easily obtain the results of Table I, where 
the location of the fixed points and the eigenvalues of the corresponding stability matrix are 
reported to leading order. 

Note that the 0(y/e) fixed points emerge at two- loop level, and are related to the de- 
generacy of the one- loop /3- functions [34]. For M = 2 the fixed points IV, VII, VIII can be 
mapped respectively into the fixed points VI, IX, X using the symmetry (2.3). 

The physically relevant fixed points are those that can be reached from the region u < 0, 
i.e. the fixed points I (Gaussian), II (O(M)), IV (Ising), VI (cubic), VII (RIM) for all values 
of M > 2, and V (mixed) and IX for M < 4 (note that in three dimensions the mixed 
fixed point V lies in the u > region for all M > 2, as shown in the preceding section). 
Concerning their stability properties, we note that all fixed points in the region u < 0, 
except the O(M) and cubic one, are unstable for any M. For M < M c = 4 — 2e + 0(e 2 ) 
(M c ~ 2.9 in three dimensions) the stable fixed point is the O(M) one, while for M > M c 
the stable fixed point is the cubic one. Note that the derivation of uoi for the fixed point VII 
(RIM) requires a three-loop calculation — we used here the results of Ref. [33]. 

In conclusion, near four dimensions the critical behavior is not changed by the addition 
of random impurities for any M > 2. Moreover, as already observed in Ref. [19] (although 
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TABLE I. Fixed points of the Hamiltonian (1.6) near four dimensions. We report the leading 
nontrivial contribution of the expansion in powers of e. Note that K d = (4vr) d r((i/2)/2. 



u/K d 



v/K d 



w/K d 



stability eigenvalues 



Gaussian 
O(M) 
SAW 
Ising 
Mixed 
Cubic 
RIM 



I 

II 
III 

IV 
V 
VI 
VII 
VIII 

IX (M > 2) 
IX, X (M = 2) 






h 



M-4 3, 
M-l 8 f 





2 v wv^ 



106 
M-4 
4(M-2) 
27 



M+8' 





2(M-1) 
M fc 



(M-2)' 



±2 A 



w u 

LO u 

w u 
|e lo u 

UJl 





2(M-4) 
3M 1 
^24 
53 1 
^24 
53" 
M-4 
3(M-2) ( 



-e, w„ 

4-M 



-e, w„ 



_ _ 4-M 

A?+8 6 ' Wl > — e ^w — M+s e 

e, u;„ = e/2, a;^ = e/2 

-e/3, = e, w w = -e/3 

„ , , 4-M , , , 4-M 
e, ^2 ~ 

M-4 
3M e 

2e, (j, 



£, UJn, 



3(M+1) 

^2 = e, ^3 



3(M+1)' 
M-4 g 
3M 



O(e), w. 

£ ' ^2 = 6 (M-2) 



^24 r 
53"V e > ^3 



53V e 



4-M 



M-4 
6(M-2)' 



that analysis missed the 0(a/c) fixed points), there is no softening of the transition for pure 
systems that are outside the attraction domain of the stable fixed point (in particular, for 
v < and for w < in the case M > M c ) . 

Even though the e expansion provides useful indications of the RG flow in lower dimen- 
sions, the validity of the extrapolation to e = 1 of the results obtained near four dimensions 
is not guaranteed, even at a qualitative level. Relevant features concerning the location and 
the stability of the fixed points can drastically change approaching three dimensions. More- 
over, new fixed points, which are not found in e-expansion analyses, may appear in three 
dimensions. For example, this occurs in the physically interesting cases of the Lp A Hamil- 
tonian describing the critical behavior of frustrated spin models with noncollinear order, 
see. e.g., Refs. [35,4], and the Ginzburg-Landau model of superconductors, where a complex 
scalar field couples to a gauge field, see, e.g., Ref. [36]. 



IV. ANALYSIS OF THE SIX-LOOP FIXED-DIMENSION PERTURB ATIVE 

EXPANSION. 

A. The six-loop series 

In the fixed- dimension field-theoretical approach, see, e.g., Ref. [4] and references therein, 
one expands in powers of appropriately defined zero-momentum quartic couplings. In the 
present case we define renormalized couplings u, v, and w from the zero-momentum four- 
point function. They are normalized so that, at tree level, they are related to the bare 
couplings u, v, and w by 

16-7T _ 167T 16n_ 

u = ——R MN um, v = ——Rmvui, w = ——wm, (4.1) 
o 00 

where R K = 9/(8 + K) and m is the renormalized mass defined by 
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rg w (p) = S ai , bj Z^ [m 2 + + 0(p 4 )] . (4.2) 

Here r^ 6 -(p) is the two-point one-particle irreducible correlation function. 

By using the numerical results compiled in Ref. [37] and a general symbolic manipulation 
program, we determined the RG functions to six loops. The resulting series for N = are 

du _ 2 2(M + 2)_ 2_ 95 , 50(M + 2)_ 2 _ 92(M + 2) 2 

Pu = m -^ — = — u + u + — — — — u v + -u w — — — u — — — — — —u v — — — — — —^uv 

' dm (M + 8) 3 216 27(M + 8) 27(M + 8) 2 

92 2 184 

„ W UW ; -' 

81 729 81(M + 8) 



50_ 2 _ 92 2 184 , _ u (u)- i - j - k \ 

u w UW ; -uvw + u > b)Auv J w , (4.3) 



\i+j+k>3 



dv _ 2 3_ 2_ 4(190 + 41 M) „ 2(131 + 25 M) 2 

= "It = — f + f H MlH « ; r-r— V : : U V 

Hv dm 2 3 27(M + 8) 2 27(M + 8) 
185_ 2 _ 400 _ 2 _ 92 2 77— / ^ ,,«■)_,■ /■ 

/ -? ) — -7 ) -7 / > — -71-7/) — ni it in _l_ -7) I \ h K i oi in 



216 u2jJ -8l(M + 8f 2w ' 729 Vw2 -8l UVW + V { J^ 3 ^ U< ^ V J ' (44) 



5v _ 2 3_ 12 308 , 104 2 832 _ 2 

p w = m- — = —w + w H — Mw + — — -vw — - — w u w —— -vw 

1 dm 2 (M + 8) 729 81 81(M + 8) 

555 2 4(370 + 23M) 2 (23M + 370) / _ , w) i , fc \ 

U w — — r-— -r-W w - — zttt r«t?w + W > of"! u l v 3 w k . (4.5) 

648 27(M + 8) 2 27(M + 8) \i + f^>3 ) 

The coefficients 6^, b^\, b\j^, with 3<i + j + A;<6, are reported in the Tables III, IV, 
and V, respectively. We have also computed the RG functions associated with the critical 
exponents to six loops. We do not report them either since they will not be used in our 
analysis (they are available on request). 



B. Resummation of the series 

The perturbative series for the M-vector model and for the cubic model are Borel 
summable and thus accurate results can be obtained by resummation methods that ex- 
ploit Borel summability and the knowledge of the large-order behavior [38,7]. On the other 
hand, perturbative series for randomly dilute models cannot be resummed naively. In the 
present case, an extension of the results of Refs. [39,40] indicates that the perturbative se- 
ries at fixed u/v and u/w are not Borel summable. For the zero-dimensional random Ising 
model, Ref. [41] showed that one could still compute correctly the free energy by means of 
a more elaborate resummation method. Although no proof exists that this procedure works 
in higher dimensions, this method was recently applied to resum six-loop expansions for the 
random Ising model in three dimensions [31], obtaining reasonably accurate results. 

Here, we apply a similar method to resum the perturbative series for the randomly dilute 
cubic model for generic values of u,v,w. If f(u,v,w) has a perturbative expansion of the 
form 
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f(u,v,w) = c n {v,w)u n , (4.6) 
c n (v,w)= J2 c nkiv k w\ (4.7) 

k,l=0 

following Ref. [31], we first resum the coefficients c n (v,w) and then, using the computed 
coefficients, we resum the series in u. The first resummation can be performed by using 
either the Pade-Borel method or the conformal-mapping method — indeed, the large-order 
behavior of the coefficients c n (v,w) is exactly the one of the M-component cubic model [7]. 
In the Pade-Borel method, for each < n < p (where p is maximum order considered), 
we consider the Pade approximants [(p — n — r n )/r n ] of the Borel-Leroy transformed series, 
which depend on an additional parameter b n . The second method uses the large-order 
behavior of the series and a conformal mapping, see Ref. [31]. Once an estimate of the 
coefficients c n (v,w) is obtained, the second resummation is performed by using the Pade- 
Borel method. We consider Pade approximants [(q — r u )/r u ] of the Borel-Leroy transformed 
series, with q < p — 1. Therefore, in the double Pade-Borel method, which we will mostly 
use, the parameters introduced by the analysis are p, {b n }, {r n } for the first resummation 
and q, b u , r u for the second one. It is impossible to vary all the b n and r n independently, 
since there are too many combinations. Therefore, following Ref. [31], we took them equal 
for all n. 

In the special cases when one of the coupling vanishes, one may use the methods already 
applied in the literature for the study of the corresponding model. For example, the series 
for u = correspond to those of the Hamiltonian (1.3); they are Borel summable, so that 
one may use the standard technique based on the knowledge of the large-order behavior 
and a conformal mapping [7]. The series for v = correspond to those of the random Ising 
Hamiltonian (1.3) with M — > 0; thus, one may use the analysis methods outlined in Ref. [31]. 



C. Results 

1. Stability of the fixed points of the pure systems 

First, we check the stability of the stable fixed point of the pure theory (u = 0) with 
respect to the perturbation induced by random dilution. 

The coordinates of the stable 0(2)-symmetric fixed point of the two-component theory 
are [27,42] u = 0, vxy = 1-402(4), and W = 0. The stability with respect to random dilution 
is controlled by 

d(3 u 



Ulu = 



du 



(4.8) 

{0,VXYfl) 



which, according to the nonpertubative argument reported in Sec. II, should be given by 

u u = (4.9) 

where axy and u X y are the critical exponents of the XY model. The analysis of the 
series (exploiting the known large-order behavior of the series and using the conformal- 
mapping method) gives tu u = 0.007(8). Therefore, the stability of the XY fixed point 
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TABLE II. Estimates of the subleading exponent uj u at the M-component cubic fixed point for 
several values of M. The last column reports the theoretical prediction —a c /u c obtained from the 
results of Ref. [7]. For M = 3 we also report the number obtained from the estimates (2.2). 



M 


[V C ,W C ] 




-OL c jv c 


3 


[1.321(18), 0.096(20)] 


0.157(28) 


0.167(24), 0.187(3) 


4 


[0.881(14), 0.639(14)] 


0.203(27) 


0.199(30) 


8 


[0.440(12), 1.136(10)] 


0.199(21) 


0.191(24) 


oo 


[0.174(6), 1.417(6)] 


0.179(16) 


0.175(30) 



is substantially confirmed, although the apparent error of the analysis does not completely 
exclude the opposite sign for oo u . The estimate of uj u is substantially consistent with Eq. (4.9). 
Indeed, olxyI v xy = —0.0217(12) (Ref. [27]), obtained from the analysis of high-temperature 
series, and [42] cyxy/vxy = —0.016(7), obtained by a more similar technique, i.e., the 
analysis of the fixed-dimension expansion of the 0(2)-symmetric model. We mention that 
perturbative studies based on shorter series [43,44] give apparently contradictory results, 
favoring a negative value of u u . 

We now consider the case M > 3. The relevant fixed point is now the cubic one with 
coordinates (0,v c ,w c ). The stability against the u perturbation is determined by the sign of 



d[3„ 



du 



(4.10) 

(0,v c ,w c ) 



Again, one expects oj u = —a c jv c , where a c and u c are the critical exponents associated 
with the cubic fixed point. Since the large-order behavior of the expansion of uj u is that of 
the series of the cubic-symmetric pure model, one can use the standard conformal-mapping 
resummation [7] . Estimates of uj u for several values of M are reported in Table II (we use the 
estimates of the cubic fixed points reported in Ref. [7]). For comparison, we also report the 
ratio —a c /u c as obtained from the results of Ref. [7]. Other estimates can be obtained using 
the results of Refs. [6,45]. For M = 3, we also quote —a c jv c fx 0.142 [45], —a c /v c fx 0.163 
[6]- 

In conclusion, the analysis of the six-loop perturbative series shows that the stable fixed 
point of the system without disorder is stable with respect to the addition of random impu- 
rities, in agreement with the Harris criterion. 



2. Fixed points induced by the disorder, i.e. for u < 0. 



We now study the stability of the RIM fixed point located in the plane v — 0. Its 
coordinates are [31] mrjm — —0.631(16), v — 0, wrim — 2.195(20). This fixed point is stable 
in the plane v — 0. To check its stability with respect to the v-perturbation we need to 
compute only 



dv 



(4.11) 



(urim,0,idrim) 
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since the derivatives of (3 V with respect to u and w vanish at v = 0. Note that the series 
for uj v is M independent. We apply the double Pade-Borel method to resum the series. We 
find that nondefective Pade approximants are obtained only for r w = 1. Setting q = p — 1, 
the results depend on four free parameters: p, r u , b u , and b w . In principle, we should look 
for the values of b u and b w that make the estimates independent of the order p of the series, 
but in the present case the results are weakly depending upon these two parameters. So, 
we consider all values < b u , b w < 10. With this choice, all approximants with r u — 1 and 
r u > 2 are defective and therefore, we use p — 4, 5, and 6 and r u = 0, 2. For r u = (direct 
summation) we obtain u v = -0.047(8){20} for p = 4, u v = 0.008(2){20} for p = 5, and 
u v = — 0.049(2){20} for p — 6, where the number between parentheses is related to the 
spread of the results of the approximants considered, while the one between braces is related 
to the uncertainty on the location of the fixed point. For r u = 2, all approximants with p = 4 
are defective. For p = 5 we obtain u v = -0.012(14){20} and for p = 6 u v = -0.080(5){20}. 
The quite large discrepancy between these estimates clearly indicates that the analysis is 
not very robust. We give as conservative final estimate uj v = —0.04(5), that includes all 
previous results [46]. This value suggests that the RIM fixed point is unstable, although it 
does not allow us to exclude the opposite case. 

Finally, we search for the presence of new fixed points in the physical region u < 0. Our 
analysis does not provide evidence for new fixed points in the physical region u < 0, at least 
for u ^ — 1, the region in which sufficiently stable results are obtained. In the case M = 2, we 
only find the fixed point predicted by the symmetry (2.3). Indeed, a fixed point equivalent 
to the RIM is expected at u* = u^ lM = -0.631(16), v* = v^ 1M + §Wr IM = 3.658(33), 
w* = — «^ IM = —2.195(20). Our resummation of the six-loop series gives consistent results 
for the location of the fixed point, i.e. u* = -0.70(5), v* = 3.7(3), w* = -2.2(2). 
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APPENDIX A: DIFFERENCES BETWEEN CUBIC AND HEISENBERG 

CRITICAL EXPONENTS 

In the three-component case the cubic critical exponents differ very little from those 
of the Heisenberg universality class. The available estimates of the critical exponents are 
indistinguishable within their uncertainty. Field-theoretical six-loop calculations for the 
cubic model [7] give v c = 0.706(6), r\ c = 0.0333(26) and 7 C = 1.390(12), while the analysis of 
Ref. [42] of six- and seven-loop series for the 0(3)-symmetric Heisenberg model provide the 
estimates v H = 0.7073(35), rj H = 0.0355(25), ^ H = 1.3895(50). We also mention the more 
accurate estimates vr = 0.7112(5), tjh = 0.0375(5), = 1.3960(9), obtained by lattice 
techniques [26]. By comparing these estimates we can only put a bound on the differences 
of the cubic and Heisenberg exponents; for example, the estimates of v differ at most by 1%. 

Much better estimates of such differences can be obtained by a more careful analysis of 
the six-loop fixed-dimension series of the cubic model (1.3) computed in Ref. [7]. We recall 
that the expansion is performed in powers of the zero-momentum quartic couplings v and 
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w associated respectively with the couplings v, w of the cubic Hamiltonian (1.3). The series 
for the cubic model can be obtained by setting u = in the series reported in Sec. IV. In 
the plane v, w the coordinates of the 0(3)-symmetric and cubic fixed points are respectively 
[42,26] X H = [1.390(4), 0] and [7] X c = [1.321(18), 0.096(20)]. In order to obtain a more 
precise determination of the exponents, one may proceed as follows. Noting the closeness of 
the fixed points, one may estimate the difference of the corresponding critical esponents r] and 
v by expanding around the Heisenberg fixed point. Therefore, the first-order approximation 
is 



Ar\ = r] c — t]h ~ A?/ 1 ) 
Av = v c - v H w Av (1) 



dv 

dv 

dv 



{v c - v H ) + 



Vh) + 



X H 



dw 

dv 
dw 



W c , 

X H 



(Al) 



X H 



The RG functions rj^v, w), r) t (v, w), and v(v, w) — (2 — 77^ + r) t ) 1 have been introduced in 
Ref. [7]. The expressions (Al) can be simplified using the relation [47] 



dw 



M + 8 drj^t 



(U,0) 



3(M + 2) dv 



(A2) 



(v,0) 



which is verified by our six-loop series and is probably exact. Therefore, we write 



Ar] 



(i) _ &n± 



dv 



Xj 



dv 



X H 



11_ 



11_ 

V C -V H + —W c 



(A3) 



In order to estimate the right-hand side of the above equations, one needs to consider 
all the elements of the covariance matrix associated with the cubic fixed point X c , i.e. 
Cuu = 0.000345, CW = 0.000380, and C^v = -0.000361, because the estimates of the 
coordinates v c ,w~c are strongly correlated. Using the covariance matrix one obtains the 
quite precise estimate v c + = 1.391(4), which is very close to vh = 1.390(4), thus 
leading to a large cancellation in the right-hand side of Eqs. (A3). In alternative, one may 
consider the approximate relation 



dP v 



dv 



(V c ~ V H ) + 



d(3 v 



X H 



dm 



0, 



(A4) 



X H 



which is obtained by expanding the equation (3 v (v c ,w c ) = 0. Within this approximation, 
one obtains v c — vh = —0.743(10) v c , and therefore v c — vh + ~^w c = —0.0009(8). Moreover, 
the analysis of the series provides the results drj ( j > /dv\ XH = 0.06(1) and dv/dv\ Xu = 0.21(1). 
Inserting in Eqs. (A3), we finally obtain A??W = -0.00005(5) and Az/W = -0.0002(2). 
Next, we determined the second-order contributions to Eq. (Al). It can be estimated by 
evaluating the second derivatives of 77^, rj t , and (3 U at X H - We obtained Ar/ 2 ) w —0.00002 
and Av^ ks —0.0001. We also checked that the third-order contributions are very small 
and negligible. Summing up, we obtain the estimates (2.1). 
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TABLE III. The coefficients b\f k , cf. Eq. (4.3). 


i i k 


M "ijk 




0,0,3 


0.090449 




0,1,2 


0.266228 + 0.00511843 M 




0,2,1 


0.23096 + 0.0403865 M 




0,3,0 


0.0513245 + 0.034637 M + 0.00448739 M 2 




1,0,2 


n ac to on 

0.467389 




1,1,1 


0.844729 + 0.0900481 M 




1,2,0 


0.281576 + 0.170804 M + 0.015008 M 2 




2,0,1 


0.857364 




2,1,0 


0.571576 + 0.285788 M 




3,0,0 


0.08992,5 




1 1 ii a 


— U.U I 0^±^±D i 




0,1,3 


—0.298245 — 0.00354189 M 




0,2,2 


—0.434652 — 0.0180279 M 




0,3,1 


-0.24561 — 0.0570548 M + 0.000878001 M 2 




0,4,0 


—0.040935 — 0.0299766 M - 0.00460823 M 2 + 0.0000731668 M 3 




l,U,o 


— U.4 / DZo4 




119 






1,2,1 


— 1.182oo — 0.250819 M + 0.004745ol M 




1,3,0 


—0.262806 — 0.187141 M - 0.0268143 M 2 + 0.000527256 M 3 




9 n 9 
z,u,z 


— L.ZZlo 




2 11 


— 2 09844 — 3441 (SfS M 




2 2 


— u.uyy4: / y — u.^u^'iui ivi — u.uo i ouuy jw 




301 


— J- . OIOOU 




310 


—0 895907 — 447Q54 M 




4A0 


-0.447316 




0,0,5 


0.0874933 




0,1,4 


0.434773 + 0.0026939 M 




0,2,3 


0.856086 + 0.0188466 M 




0,3,2 


0.805412 + 0.0689363 M + 0.000584428 M 2 




0,4,1 


0.334915 + 0.0984617 M + 0.00404184 M 2 + 0.0000476511 M 3 




u,o,u 


U.U440D04 -f- U.UoD4o0y iw -f- U.UU I lUoUZ IVI -T U.UUUZ 1 OoUy M TO.HDI4 X 1U 


6 A/f 4 
IVI 


1,0,4 


0.630396 




1 1 Q 
1,1, 


Z.404 10 + U.U0D0OOO IVI 




1,2,2 


3.52129 + 0.258766 M + 0.00231537 M 




1,3,1 


1.95691 + 0.546998 M + 0.0172536 M 2 + 0.00042021 M 3 




1,4,0 


0.326152 + 0.254242 M + 0.0484588 M 2 + 0.00150784 Af 3 + 0.0000350175 M 4 




o n q 


l.yoooy 




2 12 


5 51853 4- 296649 M 




9 9 1 

Z,Z, 1 


4.61152 + 1.18578 M + 0.0178836 M 2 




2,3,0 


1.02478 + 0.775897 M + 0.135727 M 2 + 0.00198707 M 3 




q n 9 


3.14287 




3 11 


5.26405 + 1.02168 M 




3 2 


1.75468 + 1.2179 M + 0.17028 M 2 




4A1 


2.43619 




4,1,0 


1^62413 + 0.812063 M 




5,0,0 


0.633855 




0,0,6 


-0.117951 




0,1,5 


-0.703992 - 0.00371266 M 




0,2,4 


-1.74808 - 0.0211855 M 




0,3,3 


-2.2733 - 0.0857507 M + 0.0000391428 M 2 




0,4,2 


-1.57598 - 0.192268 M - 0.00103889 M 2 + 0.0000241287 M 3 




0,5,1 


-0.517589 - 0.178763 M - 0.0115596 M 2 + 0.000202814 M 3 + 3.35858 x 10~ 6 


M 4 


0,6,0 


-0.0575099 - 0.0486174 M - 0.0112156 M 2 - 0.000619664 M 3 + 0.0000116406 M 4 + 1.86588 x 10" 7 M 5 


1,0,5 


-0.970081 




1,1,4 


-4.81696 - 0.0334495 M 




1,2,3 


-9.40257 - 0.29895 M + 0.00070522 M 2 




1,3,2 


-8.70366 - 0.994183 M - 0.00319474 M 2 + 0.00022649 M 3 




1,4,1 


-3.57656 - 1.20315 M - 0.0721365 M 2 + 0.00139263 M 3 + 0.0000411225 M 4 




1 £ n 
1,0, U 


-0.476875 - 0.398857 M - 0.0898279 M 2 - 0.00462342 M 3 + 0.000098325 M 4 


+ 2.7415 x 10~ 6 M B 


2,0,4 


-3.4864 




2,1,3 


-13.6166 - 0.329002 M 




2,2,2 


-18.9348 - 1.97645 M - 0.00712793 M 2 




2,3,1 


-10.382 - 3.38062 M - 0.185104 M 2 + 0.00215554 M 3 




2,4,0 


-1.73034 - 1.4286 M - 0.312569 M 2 - 0.0150661 M 3 + 0.000179628 M 4 




3,0,3 


-7.02194 




3,1,2 


-19.5576 - 1.50824 M 




3,2,1 


-16.0889 - 4.7868 M - 0.190108 M 2 




3,3,0 


-3.57531 - 2.85139 M - 0.574113 M 2 - 0.0211231 M 3 




4,0,2 


-8.19093 




4,1,1 


-13.4675 - 2.91439 M 




4,2,0 


-4.48915 - 3.21604 M - 0.485732 M 2 




5,0,1 


-4.86262 




5,1,0 


-3.24175 - 1.62087 M 




6,0,0 


-1.03493 
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TABLE IV. The coefficients b\f k , cf. Eq. (4.4) 



U^3 0.08O3991 + 0.0100499 M 

0,1,2 0.369295 + 0.0484367 M + 0.000284357 M 2 

0,2,1 0.602154 + 0.0995834 M + 0.00303927 M 2 

0,3,0 0.243427 + 0.0974186 M + 0.0100186 M 2 + 0.0002056 M 3 

1,0,2 0.564685 + 0.0705856 M 

1.1.1 1.80131 + 0.266837 M + 0.00520921 M 2 

1,2,0 1.0441 + 0.397331 M + 0.0362666 M 2 + 0.000364277 M 3 

2,0,1 1.44157 + 0.180197 M 

2.1.0 1.56695 + 0.524919 M + 0.0411313 M 2 
3,0,0 0.814816 + 0.101852 M 

W -0.0670637 - 0.00838297 M 

0,1,3 -0.376283 - 0.0486096 M - 0.000196771 M 2 

0,2,2 -0.857759 - 0.118296 M - 0.00138446 M 2 

0,3,1 -0.838966 - 0.178873 M - 0.00870044 M 2 + 0.0000687335 M 3 

0,4,0 -0.248229 - 0.112659 M - 0.0150567 M 2 - 0.000585531 M 3 + 2.63588 x 10~ 6 M 4 

1,0,3 -0.548393 - 0.0685491 M 

1.1.2 -2.46028 - 0.326888 M - 0.00241906 M 2 

1.2.1 -3.52974 - 0.717174 M - 0.0328982 M 2 + 0.000199551 M 3 

1.3.0 -1.35535 - 0.598611 M - 0.0763386 M 2 - 0.00279075 M 3 + 5.68096 x 10~ 6 M 4 
2,0,2 -1.86702 - 0.233378 M 

2.1.1 -5.18827 - 0.965409 M - 0.0396094 M 2 

2,2,0 -2.87884 - 1.21126 M - 0.142369 M 2 - 0.0044929 M 3 

3,0,1 -2.67974 - 0.334967 M 

3,1,0 -2.83487 - 0.985509 M - 0.0788938 M 2 

4,0,0 —1.09216 — 0.136521 M 

Tj^S 0.0777718 + 0.00972148 M 



0,1,4 0.498091 + 0.0634586 M + 0.000149661 M 2 

0,2,3 1.36139 + 0.179665 M + 0.00118632 M 2 

0,3,2 1.96208 + 0.303336 M + 0.00757647 M 2 + 0.0000396103 M 3 

0,4,1 1.35191 + 0.341883 M + 0.0225965 M 2 + 0.000151344 M 3 + 3.53109 x 10~ 6 M 4 

0,5,0 0.312662 + 0.156794 M + 0.0247779 M 2 + 0.00130674 M 3 + 6.86391 x 10~ 6 M 4 + 9.64099 x 10~ 8 M 5 

1,0,4 0.685931 + 0.0857414 M 

1.1.3 3.73862 + 0.480512 M + 0.00164801 M 2 

1,2,2 8.02393 + 1.1954 M + 0.0248927 M 2 + 0.000105193 M 3 

1.3.1 7.26959 + 1.78276 M + 0.110756 M 2 + 0.000298431 M 3 + 0.000013906 M 4 

1.4.0 2.06581 + 1.01732 M + 0.155716 M 2 + 0.00765683 M 3 + 8.67851 10" 6 M 4 + 2.53675 x 10~ 7 M 5 
2,0,3 2.68432 + 0.33554 M 

2.1.2 11.3782 + 1.62939 M + 0.0258891 M 2 

2.2.1 15.1532 + 3.58071 M + 0.208316 M 2 - 0.000313018 M 3 

2.3.0 5.60954 + 2.7006 M + 0.397804 M 2 + 0.0180874 M 3 - 0.0000496743 M 4 
3,0,2 5.6561 + 0.707013 M 

3.1.1 14.6374 + 3.06149 M + 0.153977 M 2 

3,2,0 7.87876 + 3.54347 M + 0.468803 M 2 + 0.018622 M 3 

4,0,1 5.56346 + 0.695432 M 

4.1.0 5.75768 + 2.06202 M + 0.167789 M 2 
5,0,0 1.75644 + 0.219555 M 

U^B -0.104845 -0.O131056A4 

0,1,5 -0.766484 - 0.0974606 M - 0.000206259 M 2 

0,2,4 -2.44861 - 0.314595 M - 0.00106485 M 2 

0,3,3 -4.38376 - 0.608037 M - 0.00748317 M 2 + 3.14201 x 10~ 6 M :i 

0,4,2 -4.54538 - 0.806121 M - 0.0303046 M 2 - 0.0000528742 M :i + 2.15704 x 10~ 6 M 4 

0,5,1 -2.39859 - 0.70043 M - 0.0596949 M 2 - 0.00115086 M 3 + 8.26434 x 10" 6 M 4 + 2.2795 x 10~ 7 M 5 

0,6,0 -0.453777 - 0.247702 M - 0.0450631 M 2 - 0.00309346 M 3 - 0.000053787 M 4 + 

+2.63266 x 10~ 7 M 5 + 4.89741 x 10~ 9 M 6 

1,0,5 —1.0206 — 0.127575 M 

1.1.4 -6.51597 - 0.822062 M - 0.000945689 M 2 

1.2.3 -17.4536 - 2.36625 M - 0.0229946 M 2 + 9.24619 x 10" 6 M 3 

1.3.2 -23.9918 - 4.14517 M - 0.145391 M 2 - 0.000191685 M 3 + 9.11727 x 10~ 6 M 4 

1.4.1 -15.6729 - 4.48444 M - 0.369036 M 2 - 0.00652617 M 3 + 0.000027074 M 4 + 1.11884 x 10~ 6 M 5 

1.5.0 -3.51598 - 1.89485 M - 0.337579 M 2 - 0.0223394 M 3 - 0.000357755 M 4 + 
+4.31553 x 10~ 7 M 5 + 1.52437 x 10~ 8 M 6 

2,0,4 -4.4801 - 0.560012 M 

2.1.3 -23.8866 - 3.18735 M - 0.0251901 M 2 

2.2.2 -48.8238 - 8.24755 M - 0.271114 M 2 - 0.000380338 M 3 

2.3.1 -41.9591 - 11.7782 M - 0.936933 M 2 - 0.0151248 M 3 - 0.0000114803 M 4 

2.4.0 -11.5837 - 6.16458 M - 1.07549 M 2 - 0.0685196 M 3 - 0.000997932 M 4 - 3.16906 x 10" 6 M 5 
3,0,3 -11.3503 - 1.41879 M 

3.1.2 -45.7836 - 7.2337 M - 0.188844 M 2 

3.2.1 -57.9214 - 15.4564 M - 1.13636 M 2 - 0.0136677 M 3 

3.3.0 -20.8832 - 10.7783 M - 1.79165 M 2 - 0.105048 M 3 - 0.00108945 M 4 
4,0,2 -16.8465 - 2.10582 M 

4.1.1 -41.5421 - 9.44192 M - 0.531144 M 2 

4,2,0 -21.8649 - 10.3417 M - 1.48269 M 2 - 0.066452 M 3 

5,0,1 —12.4665 — 1.55831 M 

5,1,0 -12.683 - 4.64788 M - 0.382814 M 2 

6,0,0 -3.1978 - 0.399725 M 
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TABLE V. The coefficients b\V, cf. Eq. (4.5). 



i, j, k 


M i] k 




0,0,3 


0.3t>107 




0,1,2 


1.31383 




0,2,1 


1.68533 + 0.00307141 M 




0,3,0 


0.643805 + 0.0574128 M — 0.0017162 M 




l,U,z 


1.4 / oUO 




1,1,1 


o. < yiyy + u.uuoyiuoo ivi 




1,2,0 


2.17284 + 0.193768 M — 0.00579216 M 2 




2,0,1 


2.133 




2,1,0 


o a a a a c; i ni nn77/i 7l /f 
2.44440 + 0. iyy f (4 iW 




q n n 
o,U,U 


U.yiDDDO 




4 


— (I 1 -i76 c V.>7 




n i q 


— i.ou / iy 




n o o 


— o.o4 ( ( Z + U.UU ( 041o4 ivl 




0,3,1 


—2.73858 — 0.0492189 M — 0.0000262347 M 




0,4,0 


—0.767062 — 0.0890547 M + 0.0000407114 M A — 0.0000875861 M° 




1 n ^ 


— z.uoouy 




112 


— 7 ^937 -1- fl 1 fiQfiQI M 




1,2,1 


— y. 242/2 — U.1DD114 Ivl — U.0U0Uoo0421 IVI 




1,3,0 


—3.45178 — 0.400746 M + 0.000183201 M 2 — 0.000394138 M 3 




9 fl 9 






2,1,1 


— 10.3981 — 0.212414 M 




2 2 


— ^ 82488 — 690232 M — 00248026 A/f 2 




301 


—3.89927 




310 


-4.36866 - 0.415187 M 




4,0,0 


-1.22868 




0,0,5 


0.495548 




0,1,4 


2.88579 




0,2,3 


6.8946 - 0.0230874 M 




0,3,2 


8.36453 + 0.00396206 M + 0.000213631 M 2 




0,4,1 


4.98655 + 0.175728 M - 0.00207184 M 2 - 0.0000193829 M 3 




5 


1 fiQfi^^ -t- n 1 ^791 ^Mxn nn9^ r >Szifi m 2 — n nnnnfii 471 ^ m 3 — ^ ^S7i 9 v 1 n — 6 


M 4 


lfiA 


3.24652 




113 


1 5 51 29 — 051 9466 M 




19 9 

i,z,z 


2o.2oUo H~ U.Uloo ( iy Ivl -+- U.UUU / zlUUO 7w 




1,3,1 


22.4395 + 0.790776 M — 0.00932327 — 0.0000872231 M° 




1,4,0 


6.16801 + 0.88826 M + 0.0132664 M 2 — 0.000345776 M 3 — 0.0000303026 M 4 




2 3 


8.72598 




2^2 


31.7591 + 0.124994 M 




2 2 1 


37 8666 4- 1 47045 M — 01 46658 Af 2 

1 . < !VI' IVI | J_.^|^I^|^IVJ U.lllWllllll ill 




9 q n 


lo.oio + z.U49y4 M + U.Uo4o / Vo ivl — U.UUUo49Zzl IVI 




3 2 


11.9097 




3,1,1 


28.4 + 0.965884 M 




3,2,0 


15 6128 + 2 22247 M + 0408606 M 2 




4,0,1 


7.98749 




4,1,0 


8.78218 + 0.933899 M 




5,0,0 


1.97599 




0,0,6 


-0.749689 




0,1,5 


-5.12987 




0,2,4 


— 14.928 + 0.0486813 M 




0,3,3 


-23.5697 + 0.0957169 M - 0.00083904 M 2 




0,4,2 


-21.0735 - 0.166287 M - 0.0000148277 M 2 + 4.49885 x 10~ 6 M 3 




0,5,1 


-9.82983 - 0.53385 M + 0.00220333 M 2 - 0.000130668 M 3 - 2.59594 x 10~ 6 


M 4 


0,6,0 


-1.77455 - 0.304043 M - 0.00943381 M 2 + 0.0000669939 M :i - 6.57249 x 10" 


6 M 4 - 3.75311 x 10- 7 M 5 


1,0,5 


-5.77111 




1,1,4 


-33.588 + 0.109533 M 




1,2,3 


-79.5478 + 0.323044 M - 0.00283176 M 2 




1,3,2 


-94.8309 - 0.748291 M - 0.0000667246 M 2 + 0.0000202448 M :i 




1,4,1 


-55.2928 - 3.0029 M + 0.0123937 M 2 - 0.000735009 M 3 - 0.0000146022 M 4 
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— 1 1 9782 — 2 05229 M — 0636782 M 2 4- 000452209 M 3 — 0000443643 M 4 


— 2 53335 x 10~ 6 M 5 

ij.uU'JUU /V 1U ivi 


2A4 


-18.8932 




O 1 Q 


— oy.4ylo U.UDOltw / M 




2,2,2 


-160.027 - 1.81428 M + 0.00813333 M 2 




2,3,1 


-124.409 - 7.20262 M + 0.023168 M 2 - 0.00136807 M 3 




2,4,0 


-33.6888 - 5.89631 M - 0.193047 M 2 + 0.00140166 M 3 - 0.000108167 M 4 




3,0,3 


—33.5592 




3,1,2 


-120.02 - 1.44746 M 




3,2,1 


-139.96 - 8.18899 M - 0.000367593 M 2 




3,3,0 


-50.5332 - 8.87281 M - 0.304904 M 2 + 0.00178302 M 3 




4,0,2 


—33.7557 




4,1,1 


-78^7274 - 3.5874 M 




4,2,0 


-42.6374 - 6.91385 M - 0.200772 M 2 




5,0,1 


-17.7137 




5,1,0 


-19.1868 - 2.21571 M 




6,0,0 


-3.59753 
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